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1. INTRODUCTION 

Power utilities and regulators worldwide are putting greater emphasis to find ways to reduce energy 
losses in power delivery system as it represents key indicator of an energy efficient system [1, 2]. In any 
power delivery system, 100% energy efficiency simply means that, all the electrical energy channeled into 
the network are delivered and completely consumed and accounted for by each consumer without any waste 
or loss of energy. However, this has never been possible. As energy flow in each power delivery component, 
a certain proportion of the energy are lost at each component in every level between the source to the load 
point. Therefore, energy efficient distribution system implies the minimization of energy losses that occur in 
the power delivery network. The level (percentage) of energy efficiency of power delivery network to deliver 
energy can be generally calculated as shown in (1). 


E iver MWh\|-E h 
ELoss 1%] ~ = sere > a : ° SUN (1) 
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Where: 

Fick = total energy losses (generation, transmission and distribution) 
FE aeietea = energy delivered (measured) at the generation system point 
Ee paiad = energy consumed and billed at end-customer 


Distribution energy losses can be broken down into technical losses (TL) and non-technical losses 
(NTL) [3, 4]. TL are associated with the inevitable and inherent loss of energy due to energized equipment 
and current flowing through resistive distribution components [5]. TL can be measured using energy meters 
or computed based on the network’s electrical properties, such as resistance, reactance, voltage, current, and 
power [6, 7]. Meanwhile, NTL are those energy units that are delivered for consumption but are not 
accounted or paid for, as a consequence of consumer pilferages, energy theft, faulty or incorrect billing [8]. 

Studies have shown that the average TL in distribution network worldwide ranges between 5% and 
10% of the total energy delivered, of which mostly coming from the distribution network [9, 10]. Distribution 
TL degrades the network performance and economic efficiency and increases system capacity, hence, 
resulting in billions of financial losses and increase in cost of investment [11-13]. The cost of TL and cost of 
investment in reducing TL can also significantly affect the consumer’s price or tariff of electricity. In some 
countries, as the privately-owned distribution network operation (DNO) do not directly bear these costs, 
the cost of TL is transferred as part of service cost to the customers. Consequently, regulatory bodies are 
imposing new regulations for DNO to reduce losses in their networks to a prescribed standard level or they 
will be penalized [14]. Also, distribution TL contributes to excessive greenhouse gases (GHG) emission, 
assuming it 1s entirely produced by fossil-fueled power generation plant [15, 16]. TL in distribution network 
can be reduced in many ways such as power factor correction, voltage optimization, network augmentation 
(e.g. increase conductor size), network reconfiguration (optimization), optimal sizing and location of DG and 
many more [9, 17, 18]. 

In a distribution network, majority of TL are contributed by feeders and transformers connected to 
the network. This paper concerns with TL in distribution feeders, hence, TL for transformers is not the scope 
of this paper. In this sense, the decrease of feeder TL which will improve the efficiency of network’s 
operation are constantly desired goals for utilities and regulators. 


1.1. Technical Losses In Distribution Feeders 

Electrical energy flows through multiple feeder sections to serve all connected loads. Sometimes, 
one or more line branches (or laterals) branch emanates from the main feeders. Each MV feeder and/or feeder 
section can be fairly short, on the order of a few or less than 1 km, or it can be as long as several tens of km, 
depending on the distance from the substation to the load point. Furthermore, the magnitude of load along 
connected at each load point (or load distribution) along feeder also varies at every demand period. 

TL in distribution feeder is primarily caused by “JR losses” which means that demand profiles 
containing large peaks lead to significantly more TL than flat demand profiles, even if the average power 
usage 1s the same. Theoretically, the power losses (MW) of a typical feeder in Figure | at any time t can be 
described using (2) and (3) [19]. The total TL for the entire feeder is then calculated as the sum of all TL 
contributed by all feeder section using (4). As shown in (2) is useful in calculating power (MW) losses only 
when all the operating condition and the detail parameters for each feeder section are known. 





Figure 1. One line diagram of a typical radial feeder 
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Loss,fdr _ fT Loss, fdr _ T 
E; ~~ (ey (t)dt = R; x = I? (t) (3) 
Loss,fd T 
E/°*!@" (total) = DL, (Ri x S-_, PC) (4) 
Where: 
peer af = power losses of the i feeder section 
E : peniiier = energy losses in feeders of the i" feeder section 
E : ene (total) =total feeder energy losses 
I; = current flowing of the i feeder section 
P, = real power of the i" feeder bus/load point 
V; = voltage of the i" bus 
R; = total resistance of the i feeder section 
l = feeder section/bus number/load point 
n = total number of feeder section / busses 
T = the time period of interest 


Calculating TL using the above TL formulae is difficult since they vary at every feeder section in 
the system due to variations in conductor resistance and current, caused by different length, load distribution, 
topology of the feeders [20]. It requires voluminous data collection and rigorous effort to develop, update and 
simulate/calculate the network and feeder model and many of the operating conditions occurring within the 
period are not known beforehand [21-23]. Therefore, due to these practical difficulties, distribution feeder’s 
TL are approximated using different approaches, such as the load factor approach. 


1.2. Existing Methodologies to Determine Distribution Feeder Technical Losses 

The definitive/straight forward method to assess feeder TL is to install energy meters at strategic 
locations along feeders to record the energy in and out of the individual feeder sections or any network 
component, of which would be a costly exercise [24]. Thus, to avoid this costly exercise, computational 
modelling such as load flow approach are used to perform TL assessment. 

The use of time-interval load flow simulations are widely used to accurately analyze TL in 
distribution feeders [25-27]. However, for accurate calculation of system wide distribution network, using 
these methods require in-depth knowledge and detail modelling of the distribution feeder. As the size of the 
network grows, the computation time, data and storage requirement increases exponentially. Artificial 
intelligence techniques, such as fuzzy logic [28] and clustering algorithm [29] are also applied to solve the 
problem, which requires large datasets to be trained to estimate TL. When the complete set of networks and 
energy metering data are not available, the prevailing method to estimate TL in the distribution network is to 
use the “load loss factor” method [30]. Benchmark network were used to infer TL of large distribution 
network according to their clusters [23, 24, 31]. However, since it is unlikely that any two networks and/or 
feeders exhibit the same characteristics, the benchmarking approach to infer TL of large distribution network 
might not yield acceptable results [32]. 

To address these practical difficulties, estimation of distribution feeder TL is still an aspect of great 
interest, as shown by the number of different studies performed on this subject. There is still need for 
researchers to further develop a more effective approach to estimate distribution feeder TL in the event of 
inadequate or even unavailability of accurate network and energy metering data. This paper attempt to cover 
such a knowledge gap, through developing an analytical approach to estimated TL in feeder and feeder 
sections in the absence of detail energy meter data in feeders and reducing the use of load flow simulations. 
The idea of establishing an analytical approach to estimate TL of a radial distribution feeder (G.e. with 
unidirectional power flow) is proposed. TL for each feeder sections are evaluated on a monthly basis based 
on estimation of the load profile of the load points, peak power loss generic characteristics and load loss 
factor method. Total feeder TL are estimated as the sum of all TL contributed by each feeder section. 
This paper continues and completes the preliminary research presented by the same authors in [33], 
extending feeder TL estimation model. 


2. RESEARCH METHOD 

This paper propose a new method to estimate distribution feeder TL using analytical approach. 
The feeder model is a traditional distribution feeders (without any distributed energy resources) which 
operated in radial configurations. In this model, the load is assumed balanced in all three phases and there’s 
no TL due to current in ground cables. Figure 2 shows a general radial distribution feeder with n feeder 
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sections. As the total infeed energy flows into each feeders from the source to each load point, a proportion 
of the energy are lost in each feeder sections, mainly as I7R losses. The level of TL associated with each 
feeder sections are influenced by a number of key parameters such as feeder length (/), load segment 
composition and distribution, topology, feeder loading and cable size. Each feeder sections between two (2) 
load points have different length and carry different amount of load. 


FOX. ss IN 


“energy loss 


Vv 
| : load point (11/0.4kV transformer) 


 hajetaieamas sea teeaenuaeseenn : feeder section 


v | | : source 
3(n) 
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Figure 2. Radial distribution feeder single line model with n feeder sections 


The load profile of the feeder sections depends on the amount of load characteristic measured at 
each of the load points and can be recorded using Smart Meters at all load points. However, due to economic 
reasons, not all load points are equipped with energy meters, hence, it is often times estimated or using 
sample measurement or typical load profile. Thus, in this work, analytical equations based on feeder section 
length, peak demand, peak power loss (PPL), load factor (LF) and loss factor (LsF) are formulated to 
describe the energy loss for each feeder section as well as for the entire feeder. The LF and LsF are key to 
estimate TL in distribution feeders and are one of the prevailing method used by utilities to estimate TL in 
the distribution network due to when complete set of networks and loads data are not available [34, 30]. 
This study shall adopt this approach as part of the overall proposed method. 

The load for any i*” feeder section, (P; (t)), during the time period T can be derived by aggregating 


the load at each downstream load point (P}(t)) at all t = 0...T. For n numbers of feeder sections and load 
points, the load profile for the first feeder section can be calculated as the coincident sum of all the load 
profile at load point | to n as shown in (5).The second feeder section is derived based on coincident sum 
from load profile at load point 2 to n; so on and so forth. The 30-days energy loss (in MWh) for each i*” 
feeder section (3, ) can be estimated based on its peak power loss (”9;) , loss factor (£; ) and the time period 
(1 month = 720 hours), as shown in (6) and (7). 

To further simplify the task, a set representative composite load profile, and the corresponding LF, 
are estimated for each load point to represent the type and composition of different load segment type. 
These data are obtained from local power utility based on load survey, as shown in Figure 3. 


P(t) = P(t) + PPG) +--+ PP), forallt =0..T (5) 
3, =,xXL;,x720, foralli=0..n (6) 
L=a:F,+(-a):F? ,foralli=0..n (7) 

Where: 

Di = peak power demand of i*” feeder section 

l; = length of i‘’ feeder section 

F; = load factor of i*"feeder section 

L; = loss factor of i*” feeder section 

a = loss factor coefficient 

Ap, bp, Cz, dp = polynomial coefficients of base case PPL equation for feeder section 

9; = peak power loss for the i” feeder section 

0; = PPL correction factor due to length of the i*” feeder section 
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Si = energy losses of i” feeder section 
Sp = energy losses of feeder 
Pi (t) = aggregated 24 hours load profile for the i” feeder section 
Pi(t) = representative 24 hours load profile for the i” load point 
n = total number of feeder section 
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Figure 3. Load profile curves and LF [33] 


As shown in (2), the PPL of each feeder section (§9;) depends on its PPD in MW, which varies over 
time. Traditionally, the normal method is to run load flow simulations for each feeder to find the PPL for 
each PPD, which is inefficient. To simplify analysis, this paper propose to develop a generic equation which 
enable user to estimate the PPL of feeder section at any value of PPD. To do this, a base case PPL 
characteristic equation is developed using regression analysis of load flow results.The base case feeder 
section 1s modelled in DigSILENT Powerfactory as a single feeder connected to a single end-point load, 
as shown in Figure 4. Static load flow simulations are conducted at 10%, 30%, 50%, 70%, 90% and 100% 
loading and the PPL at each PPD 1s plotted to obtain the third order polynomial regression equation. 
A sample result conducted on a sample base 11kV feeder is shown in Figure 5. This equation is then used to 
estimate the PPL of the feeder sections for any PPD for feeder of the same length and size as the base 


case feeder. 







0.014 
0.012 
a Regression equation 
3 0.01 PPL(PPD) = 5E-06(PPD)3 + 0.0015(PPD)2 - 
- 7E-06(PPD) - 2E-07 
11kV, Al XLPE, 240mm’ gee : 
Source P 5 saa 
(slack bus) Q 
4 
. < 0.004 
End-point load sa 
(2.8 MVA 0.002 
@0.95pf) , 
<—M#_—— 0 0.5 1 1.5 2 2.5 3 
l= lkm Peak power demand (MW) 


Figure 4. Single line model of base case feeder —_ Figure 5. Example regression equation of peak power loss 
section for a single 11kV feeder section 


To improve the accuracy of estimating PPL based on PPD, the PPL regression equation then need to 
be corrected to cater for different feeder length. From [33], it is shown that feeder PPL equation is linearly 
correlated with cable length. Thus, the correction factors due to length is just established by simply 
multiplying the PPL equation with the ratio of the length of the feeder of interest (/;) to the base case feeder 
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length (l,,), as shown in (8). The coefficients a,b,c and d are the PPL coefficients for the base case feeder 
section and is obtained from the graph shown in Figure 5. This equation then is used to find the 30-days TL 
for each of feeder section with different PPD and length. Finally, the total feeder TL (3,) can be calculated 


by taking the sum of TL of each of the feeder sections (3;), as shown in (9). The next section shall present 
case studies to estimate feeder losses on three (3) type of representative feeders using the proposed method. 
Figure 6 shows the flowchart of the overall methodology proposed in this paper. 


lj 
r=, x {ap? + bp? — cp; — d} (8) 


Jp = die 3; (9) 
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Figure 6. A flow chart of the proposed distribution feeder TL estimation methodology 


3. CASE STUDY 

The proposed methodology are demonstrated through several case studies based on a set of 
representative feeders (RF) which is a generic and typical feeders found based on different types of area 
served created in reference [35]. All RF, including the base case feeder comprises of 11kV, 240mm2, 3 core, 
Al XLPE cable type. While the length of the base case feeder section is chosen to be Ikm, the length of RF 
varies, depending on the location of the area served. To simplify analysis, the loads at each load points are 
assumed under balanced condition, with power factor of 0.95 and the voltage is assumed constant along 
feeder. In addition, the 30-days energy loss also assume that, the estimated load profile of the load point does 
not vary significantly for the entire loss calculation period. These generic characteristics of the RF under 
study are summarized in Table 1. 
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The 30-days TL are calculated for each feeder using the above-mentioned method. The TL 
estimation is then validated against 15-minutes time interval load flow simulation. It is important to note that, 
despite of the generalized method proposed in this paper, the RF type presented and represents the typical 
distribution feeder found in the local power utility and cannot be prescribed over the board. For feeder of 
different type, further addition and modification of the values of the parameters (e.g. PPL regression 
equation) need to be repeated. Also, in this study, the loss factor coefficient, @ , 1s assumed to be 0.25, based 
on a study performed in reference [33]. 


Table 1. Case Study Representative Feeder Characteristics 


RF Feeder PPD in ; uae NOL OAS petits Load segment Total feeder Description of 
distribution/ (and feeder 
Type MW ae eee Sections) composition length (1)inkm_ _load area served 
1 High (>4MW) Near feeder 5 100% 4] High density 
source Commercial commercial area 
Medium Nae eeaee Mixed (70% Residential 
2, (between 10 Residential, 30% 13.7 township at 
source ; 
2~3mw) Commercial) semi urban area 
4. Soweoniys, > oteenmates 10 100% residential 20.2 ae 
near feeder end villages/rural 


4. RESULTS AND DISCUSSION 

Results of the case study are organized in two (2) parts: (a) estimation of feeder TL (b) validation 
with time series load flow simulation. Table 2 and 3 shows the results of the 30 days TL estimation for each 
feeder section along the feeder as well as for the total feeder TL. TL results obtained based on the proposed 
approach shows consistency with changes in the feeder section’s characteristics. The TL of the feeder and 
feeder sections can be ranked accordingly, hence, it provides useful information to assist network planners to 
prioritize TL mitigation plan. It can be see that, although the average feeder length are short, the highest 
losses is observed in feeder section | as it has the highest PPD and LF. For RF type 2, even though the PPD 
for feeder section | and 2 are close to the rest of the feeder sections, the high losses of the two feeder sections 
are mainly contributed by high PPD and longer feeder length. The length of feeder section 10 is the highest 
among all in the same feeder but, the TL is lowest due to low PPD. From Table 3, the results of the total 
feeder TL when compared with time series load flow simulations shows less than 5% in differences. 
This shows that the proposed method yields reasonably accurate results. 


Table 2. TL Estimation for Case Studies 


cd Feeder Length , ae . ; LF LsE* Estimated branch PPL of feeder 30-days TL 
YPe section (km) ‘ane ; section PPD (MW) section (MW) (MWh) 
1 0.8 0.93 0.575 0.367 4.044 0.019515 5.50 
2 1 0.89 0.562 0.352 3.130 0.014560 3.95 
1 3 0.6 0.87 0.541 0.330 2.276 0.004604 1.17 
4 1 0.7 0.513 0.301 1.410 0.002932 0.69 
5 0.7 0.75 0.516 0.304 0.750 0.000578 0.14 
1 0.9 0.25 0.635 0.438 2.294 0.007016 233 
2 2.1 0.35 0.656 0.465 2.049 0.013041 4.57 
3 1.1 0.22 0.594 0.389 1.849 0.005562 1.65 
4 1.5 0.17 0.609 0.407 1.641 0.005969 1.85 
) 5 i 0.32 0.570 0.361 1.545 0.004227 1.17 
6 1.3 0.18 0.599 0.394 L232 0.002908 0.83 
i 0.8 0.41 0.620 0.420 1.061 0.001324 0.40 
8 1.6 0.29 0.460 0.249 0.802 0.001511 0.27 
9 Lf 0.23 0.452 0.241 0.532 0.000704 0.12 
10 iis) 0.31 0.445 0.235 0.310 0.000209 0.04 
1 12 0.11 0.465 0.254 1.808 0.005796 1.16 
2 1.4 0.13 0.466 0.254 1.700 0.005975 1.20 
3 1s) 0.11 0.464 0.252 1.583 0.004812 0.96 
4 25 0.12 0.464 0.252 1.477 0.008048 1.61 
3 5 ZA 0.11 0.464 0.253 1.361 0.005736 1.15 
6 1.9 0.1 0.466 0.254 1.251 0.004382 0.80 
w 2.6 0.25 0.466 0.254 1.154 0.005102 0.93 
8 2.6 0.32 0.461 0.250 0.925 0.003273 0.59 
9 1.9 0.33 0.453 0.242 0.629 0.001101 0.19 
10 2.7 0.31 0.445 0.235 0.310 0.000377 0.06 


Indonesian J Elec Eng & Comp Sci, Vol. 16, No. 3, December 2019: 1126-1135 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1133 


Table 3. Validation of TL Estimation with Time Series Load Flow Simulations 


RF type Estimated total feeder Simulated feeder TL Difference % 
TL (MWh) (time series) (MWh) (MWh) Difference 
1 11.45 11.87 0.42 3.51 
13.23 12.74 0.49 3.87 
3 8.66 8.94 0.29 3.19 


5. CONCLUSION 

Knowing that distribution feeders represent a significant contribution to TL, it is necessary for 
utilities to obtain a clear information on its level, location and cause of TL. A new and effective method for 
estimation TL of distribution feeder are proposed. The proposed analytical approach of using representative 
load profile and load loss factor formulation to estimate TL in each feeder sections is shown to be efficient. 
It is also robust and could be easily modified to perform TL estimation for different types and configurations 
of distribution feeders. 

The case study shows that, TL for large number of distribution feeder could be efficiently estimated 
on a monthly basis using minimal load flow simulation and minimal type of data required. TL results of 
distribution feeders estimated on a regular basis are useful to monitor the TL trend in the feeders. Operational 
plan in terms of reconfiguring the network by changing normal off points could be formulated to minimize 
TL of distribution feeders. Additionally, TL in the distribution feeders could be useful input for decision 
making in network augmentation. 

Future research work could be extended to include estimation of feeder TL with bidirectional power 
flow due to penetration of distributed generation, effect of harmonics and unbalanced condition. In addition, 
future research could also consider TL in ground cables as well as the effect of laterals to the feeder TL 
calculation. This methodology is useful in countries that Smart Meters are far from reality and the resources 
are scarce, so the system wide estimation could be performed with reasonably accurate results efficiently 
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